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ABSTRACT 

The population of gamma-ray pulsars, including Crab observed in the TeV range, 
and Vela detected above 50 GeV, challenges existing models of pulsed high-energy 
emission. Such models should be universally applicable, yet they should account for 
spectral differences among the pulsars. 

We show that the gamma-ray emission of Crab and Vela can be explained by 
synchrotron radiation from the current sheet of a striped wind, expanding with a 
modest Lorentz factor T < 100 in the Crab case, and T < 50 in the Vela case. In the 
Crab spectrum a new synchrotron self-Compton component is expected to be detected 
by the upcoming experiment CTA. 

We suggest that the gamma-ray spectrum directly probes the physics of relativistic 
magnetic reconnection in the striped wind. In the most energetic pulsars, like Crab, 
with ^ 3 g 2 /P_ 2 > 0.002 (where E is the spin down power, P is the pulsar period, and 
A" = Xj x 10* in CGS units), reconnection proceeds in the radiative cooling regime and 
results in a soft power-law distribution of cooling particles; in less powerful pulsars, 
like Vela, particle energization is limited by the current sheet size, and a hard particle 
spectrum reflects the acceleration mechanism. A strict lower limit on the number 
density of radiating particles corresponds to emission close to the light cylinder, and, 
in units of the GJ density, it is > 0.5 in the Crab wind, and > 0.05 in the Vela wind. 

Key words: acceleration of particles - magnetic fields - magnetic reconnection - 
radiation mechanisms: non-thermal - pulsars: general - gamma-rays: stars. 


1 INTRODUCTION 

It is now well established th at the pulsed spectrum of the 
Crab extends up to 400 GeV llAliu et alJl201ll : lAleksic et al.l 
120121) and very recentl y new measuremen ts of TeV photons 
have been announced dZanin et al. I [ml)- Although this is 
the most extreme examp le, several other pulsar s exhibit 
emission above 25 GeV dAckermann et all 120131 ). includ- 
ing the Vela p ulsa r with reported detection above 50 GeV 
dLeung et al.l 120141 ). The Fermi-LAT data of these objects 
have been fitted with a sub-exponential cut-off or a broken 
power-law. 

Such emission has been modelled as inverse Comp¬ 
ton scattering of magnetospheric photons by the e* 
pairs created and accelerated in the magnetosphere 
dLvutikovI 120131 ) . As an alternative, wind models have 
been proposed, which attribute the gamma-ray emis¬ 
sion to the synchrotron radiation of energized plasma ei- 
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ther very close to the light cylinder dL yubarskiil 1 19961 : 
iBai fe SpitkovskJ 20ld: Arka fe Dubusl l2013l ). or in the 
far wind zone i Petril 20121 ). Inverse Compton scattering 
of the thermal stellar photons by the instantaneously ac- 
celerated wind has been also considered in this context 
dAharonian, Bogovalov fe Khangulvanll2012l ). 


In the standard wind scenario (the “striped wind” 
model), the bulk speed of the flow is roughly constant and 
highly relativistic. The entire physics of pulsed emission is 
related to s trong beaming effects and the pre sence of cur¬ 
rent sheets dKirk. Skiaeraasen fc Gallantll200il ) - narrow re¬ 
gions in the wind, which are embedded in the col d, highly 
magn etized stripes of opposite magnetic polarity dCoronitil 
Il99d) . Radiating particles are thought to be energized dur¬ 
ing reconnection in the current layer and they are thought 
to have relativis tic thermal distrib ution that satisfies a Har¬ 
ris equilibrium dPetri fe Kirkl l2005l ) . However, recent sim¬ 
ulations of relativistic reconnection in plasmas with high 
magnetization, a 7$> 1 (where <7 is the ratio of Poynting 
flux to bulk kinetic-energy flux), suggest that the ener¬ 
gized particles are accelerated to a power-law distribution 
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llSironi fe Spitkovskvl 120141 : IWerner et al.ll20l4) . Moreover, 
it has been argued that in highly magnetized, relativistic 
systems shearing of magnetic field lines leads to formation 
of nearly force-free current layers rather than hot plasma 
sheets, and therefore a force-free equilibri um is a more ap- 
propriate model than a Harris equilibrium (IGuo et al.ll2014l . 
and references therein). 

In pulsar electrodynamics, the formation of a current 
layer can start close to the neutron star surface, along the 
last open field lines, and extend far outside the light-cylinder 
into the striped wind region. We describe a current layer as 
approximately force-free up to the light cylinder, beyond 
which reconnection sets in and proceeds over many wave¬ 
lengths in the wind as it propagates. Since reconnection 
results in energy dissipation and acceleration of particles, 
the energy density in the particle content increases and, 
together with the reconnecting field in a sheet, it main¬ 
tains the pressure balance against the fields outside the 
layer. This approach can be compared with the FIDO model 
( force-free inside - dissipative o utsi de the light cylinder) of 
iKalapotharakos, Harding fe Kazanasl (l2014h . who, however, 
attribute the energy dissipation to a changing plasma con¬ 
ductivity. They have demonstrated that the model is able 
to reproduce the variety of Fermi-LAT pulsar lightcurves. 

We focus on the spectra of gamma-ray pulsars. In our 
model the spectral shape is an imprint of the acceleration 
regime, which can be determined either by the radiative 
cooling of particles or by the size of the acceleration region 
(current sheet thickness). We show that the wind scenario 
can be distinguished from other competing models by ob¬ 
servations of its unique signature: synchrotron self-Compton 
(SSC) emission extending up to tens of TeV in the spectrum 
of the Crab pulsar. 

In section [2J we describe the basic picture of our wind 
model, including the magnetic configuration and the parti¬ 
cle distribution function. Next, as typical examples of two 
acceleration regimes we investigate the emission properties 
of the winds from the Crab pulsar and the Vela pulsar, as 
explained in section [3] Conclusions are drawn in section U 


2 THE MODEL 


2.1 Magnetic field structure 


The wind description is based on the striped wind model of 
IPetril (l2013l l. valid beyond the light cylinder. In this model 
the wind (fluid) expands with a constant relativistic Lorentz 
factor r = (1 — /3 2 ) -1 / 2 in the radial direction. In the wind 
comoving frame the electric held vanishes. We formulate the 
relevant physics in this frame and all the quantities in this 
frame are denoted with a prime, while non primed quantities 
should be understood as expressed in the lab frame. 

In the fluid comoving frame the magnitude of the mag¬ 
netic held is 


o 2 d 2 ri f sin 2 8 
3 


r 2 


+ /3 2 ^ )tanh 2 (V> s /A) 


= 5(Lc tanh 2 (?/> s /A) (1) 

where ip B = cos 8 cos x + sin 8 sin \ cos [<j> — u(t — r/pc)] de¬ 
scribes the location of the current sheet of the wind, with x 
being the angle between magnetic and rotational axes (obliq¬ 
uity), w = 2n/P being the pulsar frequency, P - the pulsar 


period, and rL = c/uj the light-cylinder radius. The current 
sheet thickness is parametrised by A A, where A = 2n/3 7 ~l is 
the striped wind wavelength and A £ [0; 0.5]. 

The ordered held exists mostly in stripes, vanishing in 
the middle of a current sheet. We assume, however, that 
the current sheet is magnetized such that there is an ad¬ 
ditional held component, generated in the course of recon¬ 
nection, which can be either shearing or random, and which 
we call “turbulent”. There are two reasons for invoking this 
component: ( 1 ) a nonvanishing synchrotron emissivity in the 
middle of the current sheet (where the ordered component 
vanishes), ( 2 ) the turbulent nature of the held seems to be 
required to explain the de crease in lin ear polarization degree 
of the Crab optical pulse JPetrill2013h . 

We assume that only a fraction ea of the available mag¬ 
netic energy is dissipated during reconnection and most 
of the held survives as a turbulent held. The simplest 
model for this turbulent component proves to be B( ur = 
&oB( nax /cosh(i/) s /A), where bo = (1 — £a) 132 . This choice en¬ 
sures that in the case of no dissipation (bo = 1 ) a turbulent 
held provides a pressure balance between a current sheet 
and the held outside. The energy dissipation and particle 
energization may be interpreted as a change in the fraction 
of plasma and held c ontrib utions to the pressure balance 
dHarrison fc Neukirchll2009l 'l. which, in the wind comoving 
frame, takes the form 


zV 2 rV2 zV 2 

-°max . . -°tu: 

- =- r P H- 

8tt 8tt y 8tt 


( 2 ) 


The pressure of ultrarelativistic particles can be approxi¬ 
mated by p ~ u'/3, where v! = me 2 f 77 /( 7^7 is the energy 
density and n'( 7 ) is the distribution function. 


2.2 Particle distribution function 

So far the current sheet in a striped wind has been modelled 
by assuming a thermal particle distribution wi th a temper - 
ature determined by the synchrotron cooling (IPetril 120121 1. 
However, simulations suggest that during reconnection par¬ 
ticles are energized to a broader power-law distribution. 
Power law indices s tn 1 are attributed to the acceleration 
at the primary X-point, whereas in dices up to s sa 3 resu lt 
from acceleration in many X-points (Ijaroschek et al.ll2004l l. 
These sp ectral features depend also on the pl asma magneti¬ 
zation a (IGuo et al .1120141 : IWerner et al.ll2014ll . such that for 
the highest values of magnetization the spectral index tends 
to one, s —> 1. Radiative cooling has not been probed in sim¬ 
ulations, but we argue that in this case a cooled distribution 
may be observed if the acceleration and cooling timescales 
are shorter than the timescale of particle injection into the 
acceleration process. 

We suggest that depending on the pulsar period P and 
spin down power E, particle acceleration at reconnection 
sites in the wind proceeds in two different regimes. 

In the first - cooling regime - the particle spectrum is 
modelled as a power law with an exponential cut off: 

n ( 7 ) = no 7 _s e _737rad (3) 

7 ( ad is the particle Lorentz factor, for which the synchrotron 
cooling timescale in the total magnetic field B[ ot = ( B 12 + 
is equal to the acceleration timescale in the recon- 
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nection electric field: 

/a \ V 2 

/ I D7T6T \ 

Trad 


3.4 x 10 5 r? /2 P^ 2 e^-2-®38 1/4 


1 + sin 2 x 
rffj- 2 + sin 2 0 


(4) 

1/4 

(5) 


where we u se Bl « (27r/P)(B/c 3 ) 1 ^ 2 (l + sin 2 x) _1,/2 

JSpitkovskvll2006 ) and the notation X = Xi x 10* in CGS 
units. We have also assumed that the reconnection electric 
field has a form E' « rB 1 , where r < 1 is a constant re¬ 
connection rate. We estimate it by comparing the recon¬ 
nection timescale t' Tec ~ T riJ(rc) w i th the wind expansion 

time fexp ~ r diss /(e d rc) 42 iannio lioH), where e d -C 1 

determines the fraction of the available magnetic energy 
that is dissipated within the distance r d i ss . Thus, the re¬ 
connection rate can be related to the dissipation distance 
r f» r 2 e d rL/r d i as (the dissipation distance in units of the 
light cylinder radius, r = r d i BB /rL, is a parameter in our 
model, and we obtain it by fitting of the computed syn¬ 
chrotron flux to the pulsar spectrum, see Sect. For usu¬ 
ally invoked reconnection rates (r ~ 0.01 — 0.2) the full dis¬ 
sipation (e d « 1) would proceed over many wavelengths i n 
the wind, possibly accelerating it rtLvubarskv fc Kirkll200llf . 
The assumption e d ^ 1 allows us to neglect the wind accel¬ 
eration. 

In the second regime, acceleration of particles is limited 
by the current sheet size, thus they escape the acceleration 
region before they re ach t h e coo ling reg ime. This case has 
been investigated by IWerner et al.l (120141 '). who observe in 
their simulations formation of a power-law particle distribu¬ 
tion with a super-exponential cut-off 

n ( 7) = no7 _s e~ 7 ^ 7sl , (6) 


where 7 b1 is the Lorentz factor, for which particle gyrora- 
dius becomes equal to a fraction </ of the system size. The 
simulations have not been initialized with a force-free con¬ 
figuration, which is of interest here, and they do not include 
radiative losses, nevertheless we expect that when the escape 
effects are dominant, a generic type of particle distribution 
© is formed, which we adopt here. In our case, however, the 
current sheet contains disordered fields, and locally the ac¬ 
celeration lengthscale is comparable to the particle confine¬ 
ment scale, which we estimate as being roughly the current 
sheet thickness A A = A2nr^/3 « A2-7rrL, thus 


7si = 


TC(A2^r L ) e B( ot 
me 2 


« 2 x 10 8 £_i A-ifi^Bgg 2 


/ rfrj~ 2 + sin 2 9 
V 1 + sin 2 X 


1/2 


(7) 

( 8 ) 


Here we have taken into account that in the lightcurves the 
width of the pulses normalized to the period of the pulsar, if 
assumed to reflect the ratio of the sheet thickne ss to the wind 
wave length, implies a parameter A ss 0.1 dPetri fc Kirkl 
l2005l l. We keep it the same for each pulsar. However, a phys¬ 
ically motivated definition of the current sheet thickness is 
based on the regime in which particle acceleration proceeds. 
It depends solely on the combination of two pulsar observ¬ 
ables: its spin down power and its period, since 7( ad and 
7 b i become comparable when E^f 2 /P -2 ~ 0.002 (where ac¬ 
cording to our fits discussed below Ti ss 2.5, fi « 3). For 


instance, for the Crab pulsar E^f 2 /P -2 ~ 3 and for the Vela 
pulsar E^g 2 /P -2 ~ 0.002. They are typical examples of two 
different scenarios. 

Synchrotron emission of an accelerated particle can be 
observed if its gyroradius is larger than the electron iner- 
tial length, a scale on which the reconnection takes plac e 
djaroschek et al.l 120041 : iTreumann &; Baumioharml [2013). 
This defines the low energy cut-off y' c in the particle dis¬ 
tribution, such that in both eq. © and eq. © 7 > 7c- In a 
general form the cut-off is given by y' c = B 1 /(SnN'mc 2 ) 1 ^ 2 , 
where the particle density is 


N' = 



(9) 


and the parameter n' Q in eq. © and eq. © is determined 
by the pressure balance, eq. ©. Calculation of 7' can be 
simplified by approximating particle distributions, eq. Q 
and eq. ©, by a power law with a sharp high energy cut-off 
at 7«d or 7 b1 , respectively. 

Using eq. © and the relation N'(r) = nNcj^cr^/Tr 2 , 
where Vgj,lc = u)Bi,/2nec is the Goldreich-Julian par¬ 
ticle number density at the light cylinder, one can con¬ 
strain the number density parameter k of radiating particles 
k « ( ecPr/B') f nfrfd'y- Note that this expression deter¬ 
mines only a lower limit on the multiplicity, because the low 
energy cut-off y' c is an estimation of the lowest energy of 
radiating particles, whose emission is observed. 


3 PULSED EMISSION 
3.1 The gamma-ray spectra 

Synchrotron flux produced by a particle distribution n l ('y) 
is given by 

~l dv £p—,If < 10) 

where 2 = (2e B CI )/(3B^ ot y 12 ) and B cr = m 2 c 3 /eh = 4.4 x 
10 13 G is the quantum critical field, 

poo 

R(Z ) = z j d£K B/a {t) « 1.78a°' 297 e _z . (11) 

e and e is the photon energy normalized to the electron rest 
mass in the lab and in the comoving frame, respectively; d 
is the distance between the observer and the pulsar, and D 
is the Doppler factor of a wind (for details see lPetri fc Kirki 
d2005l) l. Introducing V = ! 7rad an< 4 using the particle spec¬ 

trum ©, we obtain 

/ C° 2 

dye -(s+0.6)lny-y-x/y , ( 12 ) 

.ih' r ad 

where x = (2e , B cr )/(3B( ot 7( 2 d ). Asymptotic behaviour of 
this integral can be calculated by the steepest descent 
method. We find an extremum yo of the integrand from the 
condition of the vanishing of its first derivative. In the expo¬ 
nential tail 7' 7radi the solution can be found analytically 
yo ~ (2a;) 1 / 3 . Next, we expand the integrand to the second 
order around yo. Assuming y' c 7( ad , we integrate a result¬ 
ing Gaussian function. An asymptotics takes the form 

eF e syn oc 2;l-3-( S +0.6)/3 e -1.9U/ 3 ( 13 ) 
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Figure 1. Several best fits to the Cra b spectrum. Th ick gray 
l ine s hows the predicted CTA sensitivity <|de Ona-Wilhelmi et al.l 
l2012t) , black dashed line is the asymptote GS- plotted with s = 
2.2. Dashed lines show SSC components, while solid lines are the 
total (synchrotron + SSC ) spectrum. Black p oints are the data 
from Kuioer et alJ l|200lh , lAbdo et al.1 (|2010h and lAleksic et al.1 

S 


Note that t his is very similar to t he spectrum oc exp(a: 0 ' 35 ) 
obtained bv lArka fc Dublin d2013l l by a fit. 

In a similar way we calculate an asymptote of a spec¬ 
trum produced by the particle distribution (ffill . In this case 
the exponent has a form f(y) = — (s + 0.6) lny — y 2 — x/y 2 
and its derivative has a maximum at yo « a; 1 / 4 . The flux 
has the following asymptotic behaviour: 

eF e syn OC x l-3-( S +0.6)/4 e -2,V2 ( 14 ) 

The asymptotes Gsi) and da, together with fitted spectra 
(see next section) are shown in Fig.|T]and Fig. [2] respectively. 


3.2 Constraining the model parameters 


There are four free parameters of the model: the Lorentz 
factor of the wind T, the dissipation distance in units of the 
light cylinder r, the particle index s and the dissipation effi¬ 
ciency £d- We constrain them by matching the synchrotron 
flux (1101) to the phase-averaged spectra. 

To compute the fluxes in our model we integrate the 
wind emissivi ty over 3D volume (see e q. (I lOll) u sing spec¬ 
tral methods ( Press. Flannery fc Teukolskvl 19861 '). Accord¬ 
ing to this approach, the integrand is expanded in a basis 
of orthogonal polynomials and integrated term by term. As 
bases we ch oose the same special functions as proposed by 
IPetril J2013l l. i.e., Chebyshev polynomials in coordinate r, 
Chebyshev polynomials in cos#, and a Fourier series in <f>. 
This choices allow us to use a fast cosine transform to cal¬ 
culate the coefficients of expansion in the polynomial bases, 
and to apply directly the Clenshaw-Curtis quadrature rules 
for computing integrals (Press et al. 1986). 

The model spectrum is calculated at the peak of the 
lightcurve, and we obtain the phase averaged flux using 
Few = FqI, where Fq is the flux at the peak maximum, 
and I is the integral area under the lightcurve (assuming 


Figure 2. Several best fits to the synchrotron spectrum of Vela. 
For a hard particle index, the SSC component is weaker by several 
orders of magnitude and it is not shown in the plot. Black dashed 
line is the asymptote ED- plotted with s = 1 .2. Black points are 
the data from Fermi-LAT IILeung et al.ll2014l ) 


that the emission peaks are normalized to 1). This approach 
is justified by the fact that the model spectra change very 
weakly with the pulsar phase (due to a weak angle depen¬ 
dence of the Doppler factor when the emission is boosted to 
the observer frame). 

In Fig. [Qwe present several fits to the Crab data from 
Fermi-LAT and MAGIC. There is no unique solution as long 
as the measurements in the TeV range are missing. Com¬ 
puted spectra result from an interplay between all four pa¬ 
rameters of the model, which must be found simultaneously 
in order to give a reasonable fit. The green, blue and or¬ 
ange curves show the spectra with s = 2.2 and different £d, 
for which we have found the best-fit F and r. In general, 
smaller £ d lead to smaller f and larger T, thus, by increas¬ 
ing F and decreasing e d we are able to fit the data with 
the model of emission close to the pulsar, at f = 2 (orange 
curve). This fit gives an estimate of the maximum possible 
Lorentz factor of the wind, on the order of T ~ 100, similar 
to the value assumed in previous studies (lArka fc K i rldEffllH : 
iMochol &; Kirkll20i3l 'l. Note that for the best fits with the 
same value of s the brightness of the SSC component does 
not change significantly, but the maximum energy of SSC 
photons does. It is estimated by the Klein-Nishina limit, 
boosted to the observer frame: 

e « 2F 7 ; ad mc 2 w 3.6 r^P^e^E^TeV. (15) 

Red, blue and pink curves show the fits for £ d = 0.01 
and different s, for which, as before, we have also found the 
best-fit r and r. In this case, the larger is s, the smaller are 
also r and T. The brightness of the SSC component changes 
significantly, being greater for larger s. With the currently 
available VHE data the upper limit on the particle index 
can be estimated to be s ~ 2.4 (red curve). In general, soft 
indices s > 2 are required to reproduce the gamma-ray bump 
in the Crab spectrum. The low energy cut-off in the particle 
distribution is estimated to be roughly p c ~ 50. 

A fit to the Vela spectrum requires a much harder in- 
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Crab 

Vela 

Crab 

Vela 

e d 

0.01 

0.01 

Ej = 3 X 10“ ' 

e d = 10 

r 

23 

18 

82 

40 

r 

36 

12 

2 

2 

T 

0.18 

0.26 

0.001 

0.01 

K, 

7 x 10 4 

14 

0.5 

0.05 


Table 1. Reconnection rate r and the number density k of radi¬ 
ating particles obtained for two exemplary fits to the spectra. 


dex s fa 1.2 (Fig. [2j. In this case the SSC component is 
weaker by several orders of magnitude (not shown in the 
plot), because the energy density resides mainly in the high¬ 
est energy particles, which do not produce enough of low 
energy synchrotron photons that can be upscattered to very 
high energies (VHE). The low energy cut-off in the particle 
distribution is roughly 7^ ~ 2 x 10 3 . 

The number density k of the radiating particles and the 
reconnection rate r for both pulsars are compared in TablcQ] 
for two (exemplary) best-fit sets of parameters: for the effi¬ 
ciency £d = 0.01 and for the dissipation distance r = 2. In 
all cases our fits suggest that the number density of radiat¬ 
ing particles is much smaller in the Vela wind than in the 
Crab wind, bu t the reconnection rate is higher . According 
to simulations (lTakamotdl2013l : iGuo et al .11 201 dl l this result 
suggests that the Vela wind is more magnetized. 


4 CONCLUSIONS 

The spectra of gamma-ray pulsars give insight into the 
physics of relativistic reconnection in pulsar winds. Parti¬ 
cle acceleration at the reconnection sites can proceed either 
in the radiative cooling regime (in the most powerful pul¬ 
sars with Egg IP-i > 0.002), or in the size limited regime 
if the pulsar is less energetic. The synchrotron spectra be¬ 
have asymptotically like oc exp(e 1 ^ 3 ) and oc exp(e 1 ^ 2 ), re¬ 
spectively. 

In the Crab spectrum, a new SSC component at tens of 
TeV i s expected to b e observed by the upcoming experiment 
CTA (Ide Ona-Wilhelmi et ahlfioilh The MAGIC Collabo¬ 
ration (|Zanin et al]|2014h has recently obtained new mea¬ 
surements of the Crab pulsar at ~ 2 TeV. A comparison of 
our results with the spectrum in the VHE range will give a 
constraint on the particle index s of the population acceler¬ 
ated during reconnection in the wind. 

We put an upper limit on the Lorentz factor of the 
Crab wind r < 100 and of the Vela wind T < 50, which 
are obtained if the emission occurs at 1'diss = 2rL- These 
values are in good agreemen t with previous wo rks on op¬ 
tical p olarization signatures I Petri & Kirk 20051 1 and Crab 
flares jBatv. Petri fe Zenitanill2013l h The number density of 
radiating particles in the wind depends on the dissipation 
distance, but in each case our fits to the spectra imply that 
it is smaller in the Vela wind than in the Crab wind. The 
estimated reconnection rate is larger for Vela. 

Our study implies that the current sheet structure is im¬ 
printed on the pulsar gamma-ray lightcurves. Future mod¬ 
elling of the pulse profile change with energy, together with 
the phase-resolved spectroscopy, will give much more insight 
into the reconnection dynamics in pulsar winds. 
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